The alpha question became the alpha conjecture

The alpha question . Let q[1] be a prime (greater than 3, to prevent it from being 3) such that

is it then true that ((p-1)/3)!^3 = 1 (mod p )?

(Private terminology . I took to calling the other primes ' inexplicable ' for the time being. Other later examples, when I extended my searching were: 25117 - 1 = ( 2^2 ) (3) (7) (13) (23), 60919 - 1 = (2) (3) (11) (13) (71),

92401 - 1 = ( 2^4 ) (3) ( 5^2 ) (7) (11), 115827 - 1 =( 2^2 ) (3) ( 7^2 ) (197), etc etc. )

There were two first cases to examine, with unimportant order: alpha = 1 with '+1' or '-1'.

Let's say we begin, as I did, with alpha = 1 and '-1': thus we are trying to find primes { p } made up as follows:

p = 2.3*q[1]*q[2]+1 , where q[1], q[2] are primes with q[2] = 2*q[1]-1 , ( 3 < q[1] )

and are hoping that such primes satisfy:

((p-1)/3)!^3 = 1 (mod p )

Such p 's are like this:

p = 2.3*q[1]*(2*q[1]-1)+1 = 12*q[1]^2-6*q[1]+1

which is, as I observed at the time, a quadratic whose discriminant is -12

Here I show alpha = 1 and '-1' for forming some primes { p }.

I abbreviate 'alpha' to 'a' in the following to reduce the program lines in the following, in which:

> a := 1: SIGN := -1:
count:= 0: for k from 3 to 100 do
if isprime(2^a*ithprime(k)+SIGN) and isprime(2^a*3*ithprime(k)*(2^a*ithprime(k)+SIGN)+1) then count := count+1:
q1[count] := ithprime(k):
q2[count] := 2^a*q1[count]+SIGN:
p[count] := 2^a*3*q1[count]*q2[count]+1:
fi od;
array([[q[1], q[2], p], seq([q1[k], q2[k], p[k]], k=1..count)]);

matrix([[q[1], q[2], p], [7, 13, 547], [19, 37, 421...

>

At the time at which I tested the above, I had systematically tested for ((p-1)/3)! = 1 (mod p ), for all primes up to the 30,000th (note to myself: all that is saved in a Maple worksheet that I had called ' Wilson play .mws', the name of which gives an indication of my attitude up until then). The largest such prime p I had found at that stage had been 347821 (an 'inexplicable' prime, p-1 being 2^2 *3*5*11*17*31), and thus those three last new primes - 627919, 879667, and 2310019 - in the above table were first tester candidates of my alpha observation ... ( heart beating... ).

Well, I wouldn't be here talking to you if it hadn't turned out that:

> p1 := 627919: P[1,3](p1) mod p1; # didn't produce a '1'

1372

> 1372^3 mod p1; # but its cube did

1

> p2 := 879667: P[1,3](p2) mod p2;

1

> p3 := 2310019: P[1,3](p3) mod p3;

1

>

The reader may imagine how I felt seeing those 1's...

> a := 1: SIGN := -1:
count:= 0: for k from 501 to 1000 do
if isprime(2^a*ithprime(k)+SIGN) and isprime(2^a*3*ithprime(k)*(2^a*ithprime(k)+SIGN)+1) then count := count+1:
q1[count] := ithprime(k):
q2[count] := 2^a*q1[count]+SIGN:
p[count] := 2^a*3*q1[count]*q2[count]+1:
fi od;
array([[q[1], q[2], p], seq([q1[k], q2[k], p[k]], k=1..count)]);

matrix([[q[1], q[2], p], [3607, 7213, 156103747], [...

>

Here I show alpha = 1 and '+1' for forming some primes { p }:

> a := 1: SIGN := +1:
count:= 0: for k from 3 to 100 do
if isprime(2^a*ithprime(k)+SIGN) and isprime(2^a*3*ithprime(k)*(2^a*ithprime(k)+SIGN)+1) then count := count+1:
q1[count] := ithprime(k):
q2[count] := 2^a*q1[count]+SIGN:
p[count] := 2^a*3*q1[count]*q2[count]+1:
fi od;
array([[q[1], q[2], p], seq([q1[k], q2[k], p[k]], k=1..count)]);

matrix([[q[1], q[2], p], [5, 11, 331], [29, 59, 102...

>

I tested alpha up to 10 with both '-1' and '+1' choices, and checked if ((p-1)/3)!^3 = 1 (mod p ) held, not for all the produced primes { p }, but only because of computation time constraints : for example, the prime 18711862657 (see below: alpha = 7 , q[1] = 617 , q[2] = 78977 ) took 80.5 hours to test on a 2.8GHz Pentium 4 .

Here I show alpha = 7 and '+1' for forming some primes { p }:

> a := 7: SIGN := +1:
count:= 0: for k from 3 to 200 do
if isprime(2^a*ithprime(k)+SIGN) and isprime(2^a*3*ithprime(k)*(2^a*ithprime(k)+SIGN)+1) then count := count+1:
q1[count] := ithprime(k):
q2[count] := 2^a*q1[count]+SIGN:
p[count] := 2^a*3*q1[count]*q2[count]+1:
fi od;
array([[q[1], q[2], p], seq([q1[k], q2[k], p[k]], k=1..count)]);

matrix([[q[1], q[2], p], [149, 19073, 1091280769], ...

>

Of course I knew I was onto something:

The alpha conjecture. Let q[1] be a prime (greater than 3, to prevent it from being 3) such that

then ((p-1)/3)!^3 = 1 (mod p )

But the problem then became: how can one prove this ( almost certainly tru e) conjecture? That brings me to the next natural section: